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Abstract — Oscillator phase noise (PN) is one of the major 
problems that affect the performance of communication systems. 
In this paper, a direct connection between free running oscil- 
lator measurements, in terms of measured single-side band PN 
spectrum, and the optimal communication system performance, 
in terms of the resulting error vector magnitude (EVM) due to 
PN, is mathematically derived and analyzed. First, a statistical 
model of the PN, considering the effect of white and colored 
noise sources, is derived. Then, we utilize this model to derive 
the modified Bayesian Cramer-Rao bound on PN estimation, and 
use it to find an EVM bound for the system performance. Based 
on our analysis, it is found that the influence from different noise 
regions strongly depends on the communication bandwidth, i.e., 
the symbol rate. For high symbol rate communication systems, 
cumulative PN that appears near carrier is of relatively low 
importance compared to the white PN far from carrier. Our 
results also show that 1//'^ noise is more predictable compared 
to 1//^ noise and in a fair comparison it affects the performance 
less. 

Index Terms — Phase Noise, Voltage-controlled Oscillator, Col- 
ored Phase Noise, Communication System Performance, Bayesian 
Cramer-Rao Bound, Error Vector Magnitude 



I. Introduction 

OSCILLATORS are one of the main building blocks in 
communication systems. Their role is to create a stable 
reference signal for frequency and timing synchronizations. 
Unfortunately, any real oscillator suffers from phase noise 
(PN) which under certain circumstances may be the factor 
limiting system performance. 

In the last decades, plenty of research has been conducted 
on better understanding the effects of PN in communication 
systems |fTl- ||29| . The fundamental effect of PN is a random 
rotation of the received signal constellation that may result in 
detection errors ID, ifTOl . PN also destroys the orthogonality of 
the subcarriers in orthogonal frequency division multiplexing 
(OFDM) systems, and degrades the performance by producing 
intercarrier interference [|3], 0, E, [1T21 . [UTI . Moreover, the 
capacity and performance of multiple-input multiple-output 
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(MIMO) systems may be severely degraded due to PN in the 
local oscillators ifTJl . ifTSl . ||23l . Il24l . Further, performance 
of systems with high carrier frequencies e.g., E-band (60- 
80 GHz) is more severely impacted by PN than narrowband 
systems, mainly due to the poor PN performance of high- 
frequency oscillators IfTTI . ||2TI . 

To handle the effects of PN, most communication systems 
include a phase tracker, to track and remove the PN. Per- 
formance of PN estimators/trackers is investigated in ||5], ||9l , 
Ii30j . In BDl . Il32l . the performance of a PN-affected communi- 
cation system is computed in terms of error vector magnitude 
(EVM) and 1161, im, 121, 126) have considered symbol error 
probability as the performance criterion to be improved in 
the presence of PN. However, in the communication society, 
effects of PN are normally studied using quite simple models, 
e.g, the Wiener process ||22l, JH], ED, ||33l, JS]. A true 
Wiener process does not take into account colored (correlated) 
noise sources f35\ and cannot describe frequency and time- 
domain properties of PN properly ll34l . ll36l . IIJTI . This shows 
the necessity to employ more realistic PN models in study and 
design of communication systems. 

Finding the ultimate performance of PN-affected communi- 
cation systems as a function of oscillator PN measurements is 
highly valuable for designers of communication systems when 
the goal is to optimize system performance with respect to 
cost and performance constraints. From the other perspective, 
a direct relation between PN figures and system performance 
is of a great value for the oscillator designer in order to design 
the oscillator so it performs best in its target application. 

In order to evaluate the performance of PN-affected com- 
munication systems accurately, models that precisely capture 
the characteristics of non-ideal oscillators are required. PN 
modeling has been investigated extensively in the circuits and 
systems community over the past decades ll33l . ll35l . ll38l - 
[:48 |. The authors in IMI, E2, ESI have developed models for 
the PN based on frequency measurements, where the spectrum 
is divided into a set of regions with white (uncorrelated) and 
colored (correlated) noise sources. Similar models have been 
employed in l33l . P6l to derive some statistical properties of 
PN in time domain. 

Among microwave circuit designers, spectral measurements, 
e.g., single-side band (SSB) PN spectrum is the common 
figure for characterization of oscillators. Normally SSB PN 
is plotted versus offset frequency, and the performance is 
generally benchmarked at specific offset frequencies, e.g., 
100 kHz or 1 MHz O, ED, IM- In this perspective, 
oscillators with lower content of colored noise come better out 
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Fig. 1. Phase noise PSD of a typical oscillator 



in the comparison, especially when benchmarking for offset 
frequencies close to the carrier ll49l . 

In this paper, we propose a statistical PN model taking into 
account the effect of white and colored noise sources, and uti- 
lize this model to study a typical point to point communication 
system in the presence of PN. Before using the PN model, it 
is calibrated to fit SSB PN measurements of real oscillators. 
After assuring that the model describes statistical proper- 
ties of measured PN over the communication bandwidth, an 
EVM bound for the system performance is calculated. This 
is the first time that a direct connection between oscillator 
measurements, in terms of measured oscillator spectrum, and 
the optimal communication system performance, in terms of 
EVM, is mathematically derived and analyzed. Comparing 
this bound for different PN spectra gives insight into how 
real oscillators perform in a communication system as well 
as guidelines to improve the design of oscillators. 

The organization and contribution of this paper are as 
follow: In Secini we first introduce our PN model. Thereafter, 
the system model of the considered communication system 
is introduced. In Sec. HUl we find the performance of the 
PN affected communication system in terms of EVM. To do 
so, we first drive the modified Bayesian Cramer-Rao bound 
(MBCRB) on the mean square error of the PN estimation and 
identify the required PN statistics for calculation of the bound. 
In Sec.|IV]we derive the closed-from autocorrelation functions 
of the PN increments which are the required statistics for 
calculation of the system performance. Sec. |V]is dedicated to 
the numerical simulations. First, the PN samples generation for 
a given SSB phase spectrum measurement is discussed. Later, 
the generated samples are used in a Monte-Carlo simulation 
to evaluate the accuracy of the proposed EVM bound. Then, 
we study how the EVM bound is affected by different parts 
of PN. To materialize our theoretical results, the proposed 
EVM is computed for actual measurements and observations 
are analyzed. Finally, Sec. |VI| concludes the paper 

II. System Model 

In this section, we first introduce our PN model in 
continuous-time domain. Then we present the system model 
of the considered communication system. 



A. Phase Noise Model 

In time domain, the output of a sinusoidal oscillator with 
normalized amplitude can be expressed as 

V{t) = (1 + a{t)) cos (27r/ot + , (1) 

where /o is the oscillator's central frequency, a{t) is the 
amplitude noise and 4>{t) denotes the PN BTl . The ampli- 
tude noise and PN are modeled as two independent random 
processes. It is shown in ll44l that for a frequency range close 
to the carrier frequency of the oscillator, the amplitude noise 
has insignificant effect on the output signal of the oscillator 
Moreover, amplitude noise is suppressed by the amplitude 
limiting mechanism discussed in HTI . Thus, hereinafter in this 
paper, the effect of amplitude noise is neglected and the focus 
is on the study of the PN process. 

In frequency domain, PN is most often characterized in 
terms of single-side-band (SSB) PN spectrum ||33], BTl . 
defined as 

Cif) = (2) 

-TTotal 

where P(/o + /) is the oscillator power within 1 Hz band- 
width around offset frequency / from the central frequency 
/o, and i-Totai is the total power of the oscillator For an 
ideal oscillator where the whole power is concentrated at the 
central frequency, £(/) would be a Dirac delta function at 
/ = 0, while, in reality, PN results in spreading the power 
over frequencies around /q. For the further discussions, it is 
necessary to review the relation between the measured £(/) 
and power spectral density (PSD) of PN, defined as S^{f). It 
is possible to show that at high frequency offsets, i.e., far from 
the central frequency, the amount of PN is small ll33l . B3l . 
B8l . As a result, the PN PSD is well approximated with £(/) 
found from measurements at high frequency offsets ll33l . 

Cif) « S4f) for large /. (3) 

The offset frequency range where this approximation is valid 
depends on the PN performance of the studied oscillator JST] . 
It can be shown that the final system performance is not 
sensitive to low frequency events. Thus, for low frequency 
offsets, we model S^{f) in such a way that it follows the same 
slope as of higher frequency offsets. In experimental data from 
oscillators, £(/) normally follows slopes of —30 dB/dccadc 
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Fig. 2. Oscillator phase noise generation model. 

and —20 dB/decade, until a flat noise floor is reached at 
higher frequency offsets. 

According to Demir's model ll33l . oscillator PN originates 
from the white and colored noise sources inside the oscillator 
circuitry. We follow the same methodology and model PN as 
a superposition of three independent processes 



(4) 



where ^3(i) and <j>2{t) model PN with —30 and 
—20 dB/decadc slopes that originate from integration 
of flicker noise (1//) (colored noise) and white noise, 
denoted as $3(i) and $2(t), respectiveljQ. Further, (poit) 
models the flat noise floor, also known as white PN, at higher 
offset frequencies, that originates from thermal noise and 
directly results in phase perturbations. In logarithmic scale, 
the PSD of 03 (t), (f>2{t), and 4)0 (t) can be represented as 
power-law spectrums Bsl : 

Ks K2 



where K3, K2 and Kq are the PN levels that can be found 
from the measurements (see Fig. 

Due to the integration, (f>3{t) and (f>2{t) have an cumulative 
nature ll33l . BSl . PN accumulation over the time delay T can 
be modeled as the increment phase process 



$2(r)dr, (6a) 



t~T 
t 



$3(T)dr, (6b) 



t-T 



Ut.T)^<l>2it)-cj)2{t^T) 



that has been called self-referenced PN 11521 . or the differential 
PN process 153^ in the literature and it is shown that this 
process can be accurately modeled as a zero-mean Gaussian 
process (Fig. |2]i. 

B. Communication System Model 

Consider a single carrier communication system. The trans- 
mitted signal x{t) is 



N 



(7) 



'<I>3(t) and ^2{t) are frequency perturbations inside the oscillator and 
result in (j>3{t) and </>2(t) after integration. 
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Fig. 3. Communication system model with a feedforward carrier phase 
synchronizer (5]. 



where s[n] denotes the modulated symbol from constellation C 
with average symbol energy of Eg, n is the transmitted symbol 
index, p{t) is a bandlimited square-root Nyquist shaping pulse 
function with unit-energy, and T is the symbol duration ||54l. 
The continuous-time complex- valued baseband received signal 
after downconversion, affected by the oscillator PN, can be 
written as 



(8) 



where <j){t) is the oscillator PN modeled in Sec. Ill- Al and w{t) 
is zero-mean circularly symmetric complex-valued additive 
white Gaussian noise (AWGN), that models the effect of noise 
from other components of the system. The received signal (|8} 
is passed through a matched filter p*{—t) and the output is 



N 



v{t) 



X! ■*Wp(^ -nT - t)p*{-t) 
-00 „ 1 



+ / ti;(t — t)p*(— T)dr. 



dr 



(9) 



Assuming PN does not change over the symbol duration, but 
changes from one symbol to anothei^, sampling the matched 
filter output (|9]l at nT time instances results in 



y{nT) 

that with a change in notation we have 



s\n e-" 



(10) 



(11) 



where represents the PN of the n*'* received symbol in 
digital domain that is bandlimitted after the matched filter, 
and w[n] is the filtered (bandlimitted) and sampled version of 
w{t) that is a zero-mean circularly symmetric complex-valued 
AWGN with variance cr^. Note that in this work our focus is 
on oscillator phase synchronization and other synchronization 
issues, such as time synchronization, are assumed perfect. 

III. System Performance 

In this section, we find the performance of the introduced 
communication system from the PN spectrum measurements. 
Our final result is in terms of error vector magnitude (EVM), 
which is a commonly used metric for quantifying the accuracy 
of the received signal ll55l . 1561 . As shown in Fig. |3] PN 

^The discrete Wiener PN model, that is well studied in the literature is 
motivated by this assumption (e.g., 151 flol. fill , flgl . flTI . fT9l . 1221-126 1). 
We also refer the reader to recent studies of this model where the PN variations 
over the symbol peiiod has also been taken into consideration 1271 1291 . 
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is estimated at the receiver by passing the received signal 
through a PN estimator The estimated PN, denoted as 4>[n], 
is used to de-rotate the received signal before demodulation. 
The final EVM depends on the accuracy of the PN estimation. 
In the sequel, we present a bound on the performance of PN 
estimation, based on the statistics of the PN. 

A. Background: Cramer-Rao bounds 

In order to assess the estimation performance, Cramer-Rao 
bounds (CRBs) can be utilized to give a lower bound on 
mean square error (MSE) of estimation ll57l . In case of ran- 
dom parameter estimation, e.g., PN estimation, the Bayesian 
Cramer-Rao bound (BCRB) gives a tight lower bound on 
the MSE ||58l . Consider a burst-transmission system, where 
a sequence of N symbols s = [s[l], ■ ■ ■ , s[N]]'^ is transmitted 
in each burst. According to our system model (ITTl . a frame of 
signals y = . . . , y[Af]]-^ is received at the receiver with 
the phase distorted by a vector of oscillator PN denoted as 
(f = [(/)[!],..., (/)[7V]]-^, with the probability density function 
f{(p)- The BCRB satisfies the following inequality over the 
MSE of PN estimation: 



> B 



B = E^[F(cp)]+E^[-Vpog/(c^)] , (12) 

where cp denotes an estimator of (p, B is the Bayesian 
information matrix (BIM) and ">" should be interpreted as 



meaning that 



B is positive 



(13) 



semi-definite. Here, F{ip) is defined as 

F(<p) =E, [Ey|^,, [-V^log/(y|<^,s)]] 

and it is called modified Fisher information matrix (FIM) in 
the literature, and bound calculated from (fTZt is equivalently 
called the modified Bayesian Cramer-Rao bound (MBCRB) 
1591 . Based on the definition of the bound in ( fT2] i. the diagonal 
elements of B^^ bound the variance of estimation error of the 
elements of vector (p 



:^E 



(0[n] - cj)[n]f 



> 



(14) 



From (fT2]i-(fT4]i. we note that the estimation error variance 
is entirely determined by the prior probability density function 
(pdf) of the PN f{ip) and the conditional pdf of the received 
signal y given the PN and transmitted signal /(y|c^,s) (usu- 
ally denoted as the likelihood of ip). In the following, we 
derive those pdfs based on our models in Sec. Ulland use them 
in our calculations. 

B. Calculation of the bound 

1) Calculation of&^p [— log /(tp)] .• Based on our PN 
model dUi and the phase increment process defined in (|6), the 
sampled PN after the matched filter can be written as 



= 03 [n] + 02 [n] + 

n 

= 03[1]+EC3W 



00 W, 



i[l]+EC2W+0oN, (15) 



where CsW = CsinT^T) and C2N = C,2{nT,T) are the 
discrete-time phase increment processes, and 03 [1] and 02 [1] 
are the cumulative PN of the first symbol in the block, which 
are modeled as zero-mean Gaussian random variables with 
a high varianc^ denoted as c^^jj^j and c^, [1], respectively. 



According to ( fTSl l and due to the fact that (^3 [n] and C2 [n] are 
samples from zero-mean Gaussian random processes, ip has 
a zero-mean multivariate Gaussian prior /(<p) = J^iv, 0, C), 
where C denotes the covariance matrix whose elements are 
computed in Appendix |A] as 

I k 
m=2 m'=2 



from 03 [n] 
I k 

51 ^ RC2[m- m' 



m—2 m'—2 



from (/)2[?^] 



from 00 [n] 

l,k = {l...N}, 



(16) 



where i?^3 [m] and i?^2 ["^1 are the autocorrelation functions of 
Csln] and C,2[n'], and cr^^ is the variance of 0oM- The required 
statistics, i.e., R(^^[m], RQ^lm] and cr|^ can be computed 
from the oscillator PN measurements. To keep the flow of 
this section, we derive these statistics in Sec. IIVI where the 
final results are presented in (l3Ti . (l3Fi and (l4ll . Finally, 
based on the definition of f{tp), it is straightforward to show 
that V^log/(t^) = — C~^, and consequently due to the 
independence of C from (p 



[-V^log/M] -C- 



(17) 



2) Calculation of E^p [F{ip)]: According to the system 
model in (fTTl i. the likelihood function is written as 



N 



n=l 



fiy\^,s)=Y[fiy[nMnl.s[n]) 

N 

n 



^^^^{.N.*Ne-*Hj^ (18) 



where the first equality is due to independence of the AWGN 
samples. We can easily show that log /(y|t^, s) is a 
diagonal matrix where its diagonal elements are 



[V^log./(y|<P,s)]_ = 



log f{y[n]\(l)[n],s[n]) 



dcj)^ [n] 



}■ (19) 



'We consider a flat non-informative prior 1301 . 11571 for the initial PN values. 
To simplify the deiivations, it is modeled by a Gaussian distribution with a 
high variance that is wrapped to a flat prior over [0, 2it]. 
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Following (fT3] l and ( fT9] l, diagonal elements of FIM are com- 
puted as 



(20) M 



and consequently 

k=l 



\sk - sue 



1 



M 



k=l 



sin^(f) 



(27) 



where i?, is the average energy of the signal constellation. expectation in m can be computed as 

This implies that 



2Es 



-I, 



(21) 



sin^(|) 



where I is the identity matrix. Finally, from (fT2l i. ([TtI i. and 
(EB 



sin2(-)/(£)d£ 

1 



= (1 - e-'^-/2)/2, 



(28) 



B 



2^. 



where /(e) is the Gaussian pdf of e[n] as defined before. 
(22) Finally, EVM can be computed from (|25] |. (|27] | and (|28] | as 



(29) 



The minimum MSE of PN estimation ( fT4l i depends on SSB PN 
spectrum measurements through B and C^^. We will use this 

result in the following subsection to calculate a more practical As can be seen in ( |29] l, EVM is only a function of cr^[n], i.e., 
performance measure that is called EVM. 



EVM[n] = -y/2(l - e-'^|["l/2). 

1 in ( |29] l, EVM is only a fu 
the variance of the residual phase errors. 



C. Calculation of Error Vector Magnitude 

The modulation accuracy can be quantified by the EVM, 
defined as the root-mean square error between the transmitted 
and received symbols 



EVM[nl 




\Sk\n\ 



(23) 



where Sfc[n], k G {1, . . . , Af}, is the transmitted symbol from 
the constellation C with order M, at the n*'' time instance, 
and s'f.[n] is the distorted signal at the receiver. Even with 
optimal PN estimators, we have residual phase errors. Hence, 
cancellation of PN by de-rotation of the received signal with 
the estimated PN results in a distorted signal 



SkM = Sk[n\e 



i(0N-0[«]) 



(24) 



where e[n] is the residual phase error Before going further, 
assume we have used an PN estimator 1571 that reaches the 
computed MBCRB, and estimation error e[n] is a zero-mean 
Gaussian random variable. Our numerical evaluations in the 
result section support the existence of such estimators (Fig.©. 
This implies that /(e[n]) = A/'(e[n]; 0, (T^[n]), where <j1\n\ 
is defined in (fl4] | and can be computed from the derived 
MBCRB. For the sake of notational simplicity, we drop the 
time index n in the following calculations. Averaging over all 
possible values of e[n], ( |23] l is rewritten as 



EVM[n] = 



M l^k=\ 









(25) 



The magnitude square of the error vector for a given e [n] , and 
Sfc[n] is determined as 



|sfe-s;oe^"P = 2|sfe|2(l-cos(e)) 
= 4kfepsin2(|), 



(26) 



IV. Phase Noise Statistics 

As we can see in Sec. [nil the final system performance 
computed in terms of EVM ( [29] l depends on the minimum 
MSE of the PN estimation defined in ( fT4b . According to 
( fT6l l. in order to find the minimum PN variance, we have 
to compute the required PN statistics; i.e., Ri^^\m\, R(^^\m\ 
and (T? . To find these statistics we need to start from our 

00 

continuous-time PN model described in Sec. III-AI Based on 
(|6]l, 03 (t) and 4>2ii) result from integration of noise sources 
inside the oscillator On the other hand, 0o(i) has external 
sources. Therefore, we separately study the statistics of these 
two parts of the PN. 



A. Calculation of af 



00 



The PSD of 00 (^) is defined as 

5'0o(./) = Ko, 



(30) 



where Kq is the level of the noise floor that according to (|2]l, 
where we normalized the spectrum with the oscillator power, 
it is generally inversely proportional to the signal power of the 
oscillator 1381 and can be found from the measurements. We 
are interested in the part of the PN process inside the system 
bandwidth l60ll . The system bandwidth is equal to the symbol 
rate 1 /T. The variance of the bandlimited (j>o{t) is calculated 



+ 1/2T 

^0o(/)d/-^. 

-1/2T -'- 



(31) 



As 00 (0 is bandlimited, we can sample it without any aliasing. 

*This bandwidth corresponds to using a raised-cosine pulse shaping filter 
p(t) defined in (?} with zero excess bandwidth. For the general case, the 
bandwidth becomes (1 + a)/T where a denotes the excess bandwidth 1541 . 



6 



IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS-I: REGULAR PAPERS, VOL. XX, NO. X, X XXXX 



B. Calculation of i?^, [m] and i?^2 [m] 

It is possible to show that (^3 {t, T) and C2 {t, T) defined in (|6]l 
are stationary processes and their variance over the time delay 
T, is proportional to T and T^, respectively IIS, fill, 1461 . 
However, as shown in this work, their variance is not enough to 
judge the effect of using a noisy oscillator on the performance 
of a communication system, and hence their autocorrelation 
functions must be also taken into consideration. Samples of 
C3(t, r) and <^2{t,T) can be found by applying a delay- 
difference operator on (f>3{t) and (j)2{t), respectively ll33l . 1461 . 
which is a linear time invariant sampling system with impulse 
response of 

h(t.)^5{t)-5{t-T). (32) 
Starting from 02 (i), the PSD of C,2{t,T) can be computed as 
ScAI)^SUI)\H{]2^f)\\ (33) 
where H{j2Trf) = 1 — e^^'^'^f'^ is the frequency response of 



the delay-difference operator introduced in (132) . The autocor- 
relation function of (^2{t,T) can be computed by taking the 
inverse Fourier transform of its PSD 

/ + CXj 
(34) 
-00 

where r is the time lag parameter Using ( |33] l and (|34l the 
continuous-time auto correlation function can be found as 

RC2 (t) = 8 / S^, if) sm{n fTf cos(2^/T)d/. (35) 

As can be seen in ( l35T l, in order to find the closed-form 
autocorrelation functions, we do not confine our calculations 
inside the system bandwidth 1 /T. However, we see from the 
measurements that parts of (j)3{t) and 02 (^) outside bandwidth 
are almost negligible and do not have any significant effect on 
the calculated autocorrelation functions. 
The PSD of 02 (i) has the form of 

K2 



-5*02 (/) - 



(36) 



where K2 and can be found from the measurements (Fig. |3]l. 
According to ( |35] | and (|36] |. autocorrelation function of 
(^2{t,T) can be determined as 



RC2 ir)^8 -I sin(7r/T)2 cos(2^/r)d/ 

Jo J 

^ 2K2n^ {~2\t\ + \t-T\ + \t + T\) . 



Sampling (|37] | results in 



4A'27r2r if m = 







otherwise 



(37) 



(38) 



where Ri^^lm] = Rt^^ijnT). This result shows the samples 
of C2['n] are uncorrected with each other Thus, in presence 
of only 02 (i), the variance of C,2\n] is a sufficient statistic to 
describe this zero-mean process. 

To find Rc^^[m], one can consider the PSD of 03 (i) to be 
S^^{f) oc 1/ f^- However, by doing so, S(i^{f) defined in (|33l l 
diverges to infinity at zero offset frequency and hence makes it 



impossible to find the autocorrelation function in this case. To 
resolve the divergence problem, we follow a similar approach 
to l33l . l46l and introduce a low cutoff frequency 73 below 
which S^^{f) flattens. From a practical point of view, 73 is 
the lowest observable frequency and its value is proportional to 
the duration of the observation or measurement. Our numerical 
studies show that as long as 73 is chosen reasonably small, its 
value does not have any significant effect on the final result. 
Hence, we define the PSD of 03 (i) as 

K3 



(39) 



where K3 can be found from the measurements (Fig. |3]l. 
Following the same procedure of calculating i?^2(''') ( l33l - 
l35T l and using ( [39] l, the autocorrelation function of Cait) can 
be computed by solving the following integral 



K: 



P 



— sin(7r/T)2 cos(27r/r)d/. (40) 

73 



This integral is solved in the Appendix IbI Finally, the closed- 
form sampled autocorrelation function of C,z[n\ is derived as 

i?C3 [0] = -SKsn^T^ (A + log(2^73T)) (41a) 
i?C3[±l] = -SK-iTT^T^A + log(87r73T)), (41b) 

otherwise 



i?f3 [m] = - SKstt'^T'^ 
(m + 1)2 

(m - 1)2 



m2(A + log(27r73r|m|)) 
(A + log(27r73r|m+l|)) 
(A + log(27r73r|m-l|)) 



(41c) 



where A = 3/2, and F w 0.5772 is the Euler-Mascheroni's 
constant (611. The calculated variance i?(^3[0] is almost pro- 
portional to which is similar to the results of ll33l . Il46l . 
As it can be seen from i4T[ . samples of Csln] are correlated 
in this case which is in contrast to Consequently, in 

presence of 03 (t), variance of CaM is not adequate to judge 
the behavior of the oscillator in a system; it is necessary to 
incorporate the correlation properties of (sin] samples. 

V. Numerical and Measurement Results 

In this section, first the analytical results obtained in the 
previous sections are evaluated by performing Monte-Carlo 
simulations. Then, the proposed EVM bound is used to quan- 
tify the system performance for a given SSB PN measurement. 

A. Phase Noise Simulation 

To evaluate our proposed EVM bound, we first study the 
generation of time-domain samples of PN that match a given 
PN SSB measurement in the frequency domain. As shown 
in (IDi, we model PN as a summation of three independent 
noise processes 03 (t), 4>2{t), and 0o(i)- The same model is 
followed to generate time-domain samples of the total PN 
process (Fig. |4]i. Generating the samples of power-law noise 
with PSD of 1//" has been vastly studied in the literature 



M. REZA KHANZADI et al: CALCULATION OF THE PERFORMANCE OF COMMUNICATION SYSTEMS FROM MEASURED OSCILLATOR PHASE NOISE 



7 



Wo n 



W2[n\ 



W3 [n\ 
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<t>o[n] 



4>2[n] 



03 H 



Fig. 4. Phase noise simulator model. 



116211 . 11631 . One suggested approach in 11621 is to pass 
independent identically distributed (iid) samples of a discrete- 
time Gaussian noise process through a linear filter with the 
impulse response of 



Hiz) 



1 



(42) 



The discrete-time PSD of the generated noise samples can be 
computed as 



(43) 



S''{f)^alH{z)H{z-'), 



where cr^, is the variance of input iid Gaussian noise 
The scaling factor T is also needed to match S'^{f) to the 
continuous-time PSD from the measurements ||62]| . Fig. |4] 
shows the block diagram used for generating the total PN 
process. 

Generating the samples of (j)o [n] can be done by producing 
iid Gaussian noise samples, with a power to match the 
spectrum of the generated noise to the noise floor level of 
the measurements. With the noise floor level A'o and by using 
(|43T l, variance of the input iid Gaussian noise, denoted as cr^^ , 
becomes 

A"o = <,T^<, = (44) 

The second branch in Fig.|4]shows generation of the samples 
from 4'2{t) with l//^-like PSD. The liner filter with impulse 
response H2{z) — 1/(1— behaves similar to an integrator 
From ( l43T l, PSD of the generated noise can be computed as 



5l(/)=<[(l-^-^)(l-z)] 
= <[4sin2(^/rri 

„2 



(45) 



where tr^^ is the variance of input noise, the second equality 
comes from substituting jiufT for z, and approximation is 
valid in a frequency range lower than sampling frequency. 
Suppose /\2 is the level of 1/ P noise identified from the 
SSB PN measurements. Exploiting ( l45T l, the variance of the 
iid input noise to the first branch can be computed as 



K2 

P 



Aan^pT^ 



(46) 
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Fig. 5. PSD of generated PN samples vs. theoretical PSD. Generated phase 
PSD is matched to the desired PSD with the given values of A'o, and 



A linear filter with impulse response of H^lz) = 1/(1 — 
^-1^3/2 jg ^ggj jjjg third branch to generate the time- 
samples of 1//^ noise. Employing (l43T l and taking the same 
approach as ( |45] |. the PSD of the output can be found as 



^1(/) 



2 



[(4sin2(^/T))^ 



zn 
1-1 



« — (47) 

Finally, using (07]), the variance of iid input noise to the 
third branch is computed as 

El 

P 



tT - 



2„3 



(48) 



(27r/r)3- 

Implementation of H^{z) can be a matter of discussion. For 
more information, the reader is directed to fl62l that covers 
different aspects of implementation of such a filter 

Fig. |5] shows the total one-sided PSD of the generated PN 
samples for a particular example. The frequency figures of 
merits are set to be = -110 dB, A'3 = 10"', and K2 = 10. 
Employing (l44l i. ( |46] |. and (08]), the variance of input white 
Gaussian noises to the PN generation system in Fig. 0] for a 
system with symbol rate 10^ symbol/sec are calculated to be 
cr2 ^ = 5 X 10-^ al^ = 1.97 x 10"^ and cr^ ^ = 1.26 x 10"^ 
This figure shows that generated time-domain samples match 
to the PSD of PN. It can also be seen that the PSD of the 
simulated PN flattens below certain offset frequencies and it 
has a Gaussian shape. We refer the reader to ll43l . ll44l . ll63l for 
more details on mathematical explanation of this phenomenon. 

B. Monte-Carlo Simulation 

Consider a communication system like that of Fig. [3] 16- 
QAM modulation is used and length of the communication 
block is set to 200 symbols. A local oscillator with the PN 
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Fig. 6. The phase error distribution of the symbol in middle of block (n = 
100) compared to a Gaussian fit. It can be seen that the phase en'or distribution 
is zero-mean Gaussian with a good approximation. 
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Fig. 7. Proposed theoretical EVM bound vs. the EVM from the Monte- 
Carlo simulation over a block of 200 symbols. The PSD in Fig. [5] is 
considered as the PN PSD. 16-QAM modulation is used, symbol rate is set 
to 10^ [Symbol/sec] and variance of AWGN cr^ = IQ-^. 



PSD of Fig. [5] is used. For the Monte-Carlo simulation, we 
first generate the PN samples following the routine proposed 
in subsection I V-AI Then, we design the maximum a posteriori 
(MAP) estimator of the PN vector ip at the receiver The MAP 
estimator finds (p that maximizes the posteriori distribution of 

Vmap = argmax f{(p\y,s) 
v 

= argmax fiy\(p, s)f{(p). (49) 

The needed likelihood and prior functions for designing this 
estimator are calculated in Sec. However, the detailed 
implementation of this estimator is not in the focus of this 
paper and we focus only on the final results. The estimated 
phase values are used to eliminate the effect of PN by de- 
rotation of the received signals. Finally, the EVM is computed 
by comparing the transmitted symbols with the signal after PN 
compensation. Fig. |6] shows the density of the residual phase 
errors for one of the symbols in the frame. It can be seen 
that phase errors are almost zero mean and have Gaussian 
distribution. Fig. [T] compares the proposed theoretical EVM 
bound against the resulted EVM calculated from the Monte- 
Carlo simulation. It can be seen that the calculated EVM 
from the empirical simulation perfectly matches the proposed 
theoretical bound. 

C. Analysis of the Results 

Now, when the EVM bound is evaluated, the system perfor- 
mance for a given oscillator spectrum may be quantified. In 
this section we study how the EVM is affected by white PN 
(PN floor) and cumulative PN, respectively. The effect from 
cumulative PN is further divided into origins from white and 
colored noise sources, i.e., SSB PN slopes of —30 dB/decade 
and —20 dB/decade, respectively. It is found that the in- 
fluence from the different noise regions strongly depends on 




Noise Floor Level [dBm] 



Fig. 8. Proposed theoretical EVM bound against different noise floor levels. 
K2 = 1 and K:i = 10^ are kept constant. The low cut-off frequency 73 is 
considered to be 1 Hz and cr^ = 10^'^. 16-QAM modulation is used and 
block-length is set to 10. In the hatched regime, the white PN (noise floor) 
dominates over the cumulative part of the PN. 

the communication bandwidth, i.e., the symbol rate. For high 
symbol rates, white PN is more important compared to the 
cumulative PN that appears near carrier 

Fig. [8] compares the performance sensitivity of two com- 
munication systems with different bandwidths, namely GSM 
and W-CDMA against a set of different noise floor levels. 
GSM operates with the symbol rate of 0.271 MSymbols/s 
that leads to 3.69 /is symbol duration. In contrast, W-CDMA 
has 3.84 MSymbols/s symbol rate results in 0.26 /is symbol 
time that is almost 14 times shorter than that of GSM ll64l . 
It is seen in Fig. [8] that an increase in the level of white PN 
affects the W-CDMA system with high symbol rate much more 
than the more narrowband GSM system. This result can be 
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Fig. 9. Proposed theoretical EVM bound against different values of comer 
frequency /corner for GSM and W-CDMA systems. K2 =0.1 and Ko = 
— 160 dBm are kept constant and /corner is increased by adding to K^. The 
low cut-off frequency 73 is considered to be 1 Hz and crj, = lO^**. 16-QAM 
modulation is used and block-length is set to 10. 
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Fig. 10. Two SSB PN spectrams with pure 1//- PN and 1//^ PN. We 
assume the two spectrums have a very low white PN level. For a W-CDMA 
system, both spectrums result in the same variance of phase increments 
i?^(T = 0) = 6.2 X 10-^ [rad^]. 



intuitively understood, since in a system with a higher symbol 
rate, symbols are transmitted over a shorter period of time 
and thus experience smaller amount of cumulative PN. On the 
other hand, the amount of phase perturbation introduced by 
the white PN is a function of the system bandwidth and a 
wideband system integrates a larger amount of white PN (ISTi . 
Therefore, in contrast to the cumulative PN, white PN affects 
a system with high bandwidth more compared to a system 
with a narrower bandwidth. 

The next step is to identify the different effects from 
cumulative PN originating in white noise sources (slope 
—20 dB/decade) and cumulative PN originating in colored 
noise sources (slope —30 dB/dccadc). Fig.|9]shows the effect 
of changing the corner frequency on the performance of GSM 
and W-CDMA systems by increasing the level of 1//'^ noise, 
K^. Other parameters such as A'2 and Kq are kept constant 
in this simulation to just capture the effect of different values 
of if 3. Intuitively the performance degrades when the noise 
level is increased. However, as seen in Fig. |9] the EVM 
is not significantly affected below certain comer frequencies 
(./corner < 1 kHz for GSM and /comer < 100 kHz for W- 
CDMA). This constant EVM is due to the dominant effect of 
1/ P on the performance. By increasing the corner frequency, 
after a certain point 1/P becomes more dominant which 
results in a continuous increase in EVM. It can also be seen 
that the GSM system is more sensitive to increase of the 
1/ P noise level. Because of the higher bandwidth, the W- 
CDMA contains more of the 1/P noise which is constant 
and dominates the 1/ P effect, and its EVM stays unchanged 
for a larger range of comer frequencies. 

Finally, we compare the individual effect of 1/ P PN and 
1 //"^ PN on the performance. Consider two SSB PN spectrums 
as illustrated in Fig. [TO] One of the spectrums contains pure 



1/ P PN while 1/ P PN is dominant in another In a W- 
CDMA system with the symbol rate of 3.84 MSymbols/s 
(bandwidth of 3.84 MHz), the variance of phase increment 
process for both spectmms is equal to i?^ (r = 0) = 6.2 x 
10^^ [rad^]. However, comparing the EVM values shows that 
the spectrum with pure 1/ p PN results in 2.36 dB lower 
EVM. This is due to the correlated samples of phase increment 
process for \/ p noise which results in lower PN estimation 
errors compared to 1/ P noise. 

D. Measurements 

To materialize the analytical discussion above. Fig. [TTIand 
Fig. [12] show the measured SSB PN spectmms from a GaN 
HEMT MMIC oscillator under two different bias conditions 
with drastically different characteristics for the cumulative PN. 
Fig. [TT] shows the spectrum for the oscillator biased at a drain 
voltage of Vdd = 6 V and drain current of Id = 30 mA. 
At this bias condition, the 1// noise (flicker noise) from 
the transistor is fairly low. The comer frequency between 
the —30 dB/decade and —20 dB/decade regions can be 
clearly detected at 83.3 kHz. In contrast. Fig. [12] shows a 
spectrum from the same oscillator biased at Vdd = 30 V 
and Id = 180 mA. Under this bias condition the noise 
from colored noise sources is increased significantly and the 
cumulative PN has —30 dB/decade slope until it reaches the 
white PN floor Further, the power of the oscillator is higher 
in Fig [TT] resulting in a lower level for the white PN. Fig. [T3] 
compares EVM for the SSB PN spectrums in Figs. [TTIand [TZi 
respectively, versus the symbol rate. As expected based on the 
results in subsection IV-CI the spectmm in Fig [TT] gives the 
best EVM for low symbol rates, while the spectrum in Fig [12] 
gives the best EVM for higher symbol rates as a result of the 
lower level of white PN. 
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Fig. 11. SSB PN spectrams from a GaN HEMT MMIC oscillator. Drain 
voltage Vdd = 6 V and drain current Id = 30 mA. The comer frequency 
at /corner =83.3 kHz. 
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Fig. 12. SSB PN spectrams from a GaN HEMT MMIC oscillator. Drain 
voltage Vdd = 30 V and drain cun'ent Id = 180 mA. 



VI. Conclusions 

In this paper, a direct connection between oscillator mea- 
surements, in terms of measured single-side band PN spec- 
trum, and the optimal communication system performance, in 
terms of EVM, is mathematically derived and analyzed. First, 
we found the statistical model of the PN which considers the 
effect of white and colored noise sources. Then, we utilized 
this model to derive the modified Bayesian Cramer-Rao bound 
on PN estimation that is used to find an EVM bound for the 
system performance. 

The paper demonstrates that for high symbol rate communi- 
cation systems, the near carrier cumulative PN is of relatively 
low importance compared to white PN far from carrier Our 
results also show that 1//^ noise is more predictable compared 
to 1//^ noise, and in a fair comparison it affects the system 
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Fig. 13. EVM comparison of given measurements in Fig. II 1 1 and 1121 vs. 
symbol rate (bandwidth). The low cut-off frequency 73 is considered to be 
1 Hz, 16-QAM modulation is used and crj, « 0. 



performance less. These findings will have important effects 
on design of hardware for frequency generation as well as the 
requirements on voltage controlled oscillator design, choice of 
reference oscillators and loop bandwidth in the phase locked 
loops. 

Appendix A 

In this appendix, elements of the covariance matrix C are 
calculated. According to (T5[ 

[C];,fc=E[(#]-E[#]])(0[fc]~E 

E [ (00 [I] + 03 [1] + 02 [1] + 5] (C3 [m] + C2 H ) 

m=2 

k 

X (00 M + 03 [1] + 02 [1] + ('^sK] + C2K]) 

m'=2 

03[1]03[1]] +e[02[1]02[1]] +E[0oH0o[fc]" 

; k 

+ E E ie[c3HC3K]_ 

m=2 m'=2 

; k 

+ Rq^ [m — m'] + [m — m'] 



E 



m— 2 m'—2 

l,k = {l...N}. (50) 



Appendix B 

Steps taken to solve the integral in ( i40] i are described here. 
We can write ( i40] i as 



P + ^l 



siii(7r/T)2cos(27r/T)d/ 
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+00 



= 4 
- 2 



P + ll 
+00 J. 



cos(27r/T)d/ 



3Cos(2^/(|r + T|)d/ 

7+73 



- 2 



+ CXD 



P + 



cos(2^/(|r-T|)d/. (51) 



It is clear that solving the integral in the form of 
/o^°° f3^+^--i cos(27r/T)d/ is enough to compute the total inte- 



gral of (Dil i. This integral is complicated enough that powerful 
software such as Mathematica are not able to converge to the 
final answer Consequently, first, partial-fraction decomposi- 



tion of ,^ i is done: 



B 



C 



,-j27r/3 



pj27r/3 



(52) 



Note that, 73 is a real positive number Using Mathematica 
(Version 7.0), the following integral can be evaluated 

n+oo 2 



f + /3 



cos(27r/r)d/ = 

— cos(2/37rr)cosint(— 2/37r|r|) 

- I sin(2/37r|T|)(7r + 2 sinint(2/37r|r|)), (53) 



where f3 must be a complex or a negative real number, and 
sinint(-) and cosint(-) are sine and cosine integrals define as 

. . . N r sin(t) , 
sinint(x) = / — ■ — at 



cosint(.T) = r + log(2:) + / ^^^^^ -dt, (54) 



/o t 

where F « 0.5772 is the Euler-Mascheroni's constant. 

Consider the case where time lag r is small. By Taylor 
expansion of the functions in ( l53T l around zero 

sin(a;) ~ x — + . . . cos(.t) = 1 — + . . . 

6 2 

sinint(a;) = x V ■ ■ ■ cosint(a;) = F + log(.T) V ■ ■ 

18 4 

and neglecting the terms after second order, the integral can 
be approximated as 

n+00 -j^ 

cos(27r/r)d/ « 

- F - log(-2/37r|T|) - 2/37r^|r| 
, (2/?7r|r|)2 



/ + /? 



-(A + log(-2/3^|r|)), (55) 

where A = F — |. Employing this approximation and the 
fraction decomposition in ( |52] |. followed by a series of sim- 
plifications 

-——3 cos(27r/r)d/ sa 
Z"* + 73 



37I VV3 



672^V2(A + log(2737r|r|)) ) . (56) 



Now the first term in dSll is calculated. By changing the 
variable t to r + T and r — T, second and third terms can 
also be computed, respectively. Finally, R^z is approximated 
by 



2 

{r-Tf 



-r2(A + log(2^73|r|)) 
■(A + log(27r73|r + r|)) 
■(A + log(27r73|T-r|)) 



(57) 



To calculate the ACF for r 0, and r = |T|, we need to take 
the limits of (ISTT i as r approaches 0, and |T|, respectively that 
results in 



lim Rc^^{t) = -8A'37r2r2(A + log(27r73r)), 



(58) 



lim i?o,(T) = -8^5:3 TT^T^ (A + log(87r73r)). (59) 
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